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Abstract 



We present the results of deriving the Israel-Stewart equations of relativistic dissipative fluid 
dynamics from kinetic theory via Grad's 14-moment expansion. Working consistently to second 



^ I order in the Knudsen number, these equations contain several new terms which are absent in 



previous treatments. 



^ ' 1 Introduction 

I Fluid dynamics has always been an important qualitative tool to describe the collective flow of hot 
and dense matter created in heavy- ion collisions However, the great success of the fluid dynamical 
^ ■ model only came with nuclear collisions at RHIC energies, where for the first time fluid dynamics was 
00 able to describe flow observables, such as the elliptic flow [2], on a quantitative level [3]. 

] The respective calculations [3] of the elliptic flow at RHIC energies were based on ideal fluid dy- 
I namics, i.e., all dissipative effects were neglected. The good agreement between data and the ideal fluid 
^ ' dynamical calculations gave rise to the notion that "RHIC scientists serve up the perfect fluid". This 
■ implies that dissipative effects have to be small in order to not spoil the agreement with data. 

However, the fluid dynamical equations of motion are partial differential equations which require 
initial conditions in order to solve them uniquely. In principle, these introduce inflnitely many degrees 
of freedom that can be tuned to achieve agreement with the data (once an appropriate "freeze-out" 
prescription is adopted to convert the fluid into particles [1]). Since the initial conditions in heavy- 
ion collisions are not known sufficiently precisely, one has to perform calculations within relativistic 
dissipative fluid dynamics and for various realistic initial conditions in order to conflrm the smallness 
of dissipative effects. 

Unfortunately, a consistent, stable, and causal formulation of relativistic dissipative fluid dynamics is 
far from trivial. The so-called first-order theories due to Eckart ^ and Landau have been shown |7] 
to lead to unstable solutions and to support acausal propagation of perturbations. A viable candidate 
for a relativistic formulation of dissipative fluid dynamics, which does not have these problems (for a 
large class of equations of state [Ej), is the so-called second-order theory due to Israel and Stewart [9]. 
In recent years, theoreticians have begun to apply this theory for the description of collective flow in 
heavy- ion collisions [TOj [11], [12] . 
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In this paper we present the results of an analysis (the details of which will be presented elsewhere 
|13j ) of deriving the Israel-Stewart equations of relativistic dissipative fluid dynamics from kinetic 
theory, using the Boltzmann equation and Grad's 14-moment [Hj. We show that, working consistently 
to second order in the Knudsen number, additional terms arise in these equations that have been missed 
in the original paper [9] and in subsequent treatments flOi [TT| [T2]. 

This paper is organized as follows. After introducing our notation in Sec. [21 we discuss the power 
counting scheme in terms of the Knudsen number (Sec. [3]). In Sec. |4]we present and discuss the Israel- 
Stewart equations. A concluding section summarizes our results and gives an outlook to future work. 
Our units are h = c = ks = 1, the metric tensor is g^"^ = diag(+, — , — , — ). The scalar product of 
4-vectors A^, 5^ is denoted as A'^g^^yB'' = Af'B^ = A - B. 

2 Preliminaries 

The quantities that are evolved through the fluid dynamical equations are the 4-current of net charge, 
iV^ (for the sake of simplicity, we only consider a single conserved charge, the generalization to several 
charges is straightforward [IS]), and the symmetric rank-2 energy- momentum tensor, T^*^. In the 
following, we first discuss the tensor decomposition of these quantities and then give the corresponding 
fluid dynamical equations which express the conservation of net charge and energy-momentum. Finally, 
the Navier-Stokes approximation to relativistic dissipative fluid dynamics is discussed. 

2.1 Tensor Decomposition of andT^^ 

The tensor decomposition of the net charge 4-current with respective to an arbitrary 4-vector reads 

N^' = nu^' + . (1) 

We shall identify with the fluid ^-velocity which is time-like and normalized to one, u ■ u = 1. This 
yields = 7(1, v), where v is the fluid 3-velocity and 7 = (1— v^)~^/^ the corresponding Lorentz gamma 
factor. Due to the normalization condition, the 4-vector has only three independent components. 
The component of A^'^ in the direction of is the net charge density in the fluid rest frame, n = N- u. 
The component orthogonal to is the diffusion current, i.e., the flow of net charge relative to u^, 
z/^ = A'^'^Niy. Here, A^*^ = g'^'^ — u^u^ denotes the projector onto the 3-space orthogonal to u^. By 
construction, /S.^^Uy = 0, and thus also u -u = 0. Therefore, the 4-vector has only three independent 
components. 

The tensor decomposition of the energy-momentum tensor reads 

T^"' = e u^'u'' -{p + U) A^" + 2 q^^u"'' + n^"" , (2) 

where e = u^T^^u^ is the energy density in the fluid rest frame. The projection p + 11 = — | A^^^T^^ 
is the sum of thermodynamic pressure, p, and bulk viscous pressure, U. The heat flux current, i.e., 
the flow of energy relative to u'^, is q'^ = A^^Tyxu'^- By construction, g ■ u = 0, and q^ has only 
three independent components. The notation a'-"^ stands for symmetrization in all Lorentz in- 
dices, e.g., a'^^"'' = ^{a^^" + a"^). Finally, tx^^" = T^^''^ is the shear stress tensor, where a^^"^ = 

[A^^'^A"^!^ — I A^^^ Aapj a"^ denotes the symmetrized, traceless spatial projection. Thus, by construc- 
tion, TT^'^Ufj, = Tx^^Uy = 7r^^ = 0. This implies that tt^'^ has only five independent components. 

2.2 Fluid Dynamical Equations of Motion 
The conservation of net charge reads 

d- N = h + n9 + d = , (3) 
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where a = u - da denotes the comoving derivative. In the fluid rest frame, Uj^p = (1, 0, 0, 0), it is simply 
the time derivative, orf = dta. The quantity 6 = d ■ u is the so-called expansion scalar. 

The equation for energy-momentum conservation, d^JI"^^ = is a 4-vector and can also be tensor- 
decomposed into a component parallel to m^, 

u.d^r^^" = e + {t + p + Y])d + d-q-q-u-Ti^'' d^u, = Q , (4) 

which represents the conservation of energy, and into the three independent components orthogonal to 
u^, A^*^ d^Ty\ = 0, which can be cast into the form 

(e + p) = v^(p + n) - n m'^ - A'^^g^ - q^" 6 - q ■ du^ - A^"" d^n^x , (5) 

where the notation = A^'^d^ stands for the 3-gradient, i.e., the spatial gradient in the fluid rest 
frame. Equation ([5]) is the so-called acceleration equation. The first term on the r.h.s. expresses the 
fact that changes in the fluid 4-velocity are driven by pressure gradients. 

The problem with the five equations of motion ([3]), (jl]), and ([5]) is that, for any given fluid 4-velocity 
u'^, they contain 15 unknowns, e, p, n, 11, the three components of u'^, the three components of g'^, and 
the five components of tt^'^. The choice of reference frame does not alleviate this problem: in the Eckart 
frame z/^ = 0, and in the Landau frame q^ = 0, however, then u'^ is no longer fixed, but becomes a 
dynamical quantity, and one still has 15 unknowns. In principle, there are two strategies to continue. 

1. One assumes the fluid to be in local thermodynamical equilibrium. Then, all dissipative quantities 
vanish, 11 = (or z/'') = tt^'^ = 0, which eliminates nine out of the 15 unknowns. This is the so- 
called ideal-fluid limit. In this case, the reference frame is uniquely defined. The sixth remaining 
unknown, for instance p, is determined by choosing an equation of state for the fluid in the form 
p = p{e,n). Then, the system of five equations ([3]), (jl]), and ([5]) can be uniquely solved for e, n, 
and u'^. 

2. One provides equations that determine the dissipative quantities: the equations of dissipative fluid 
dynamics. One distinguishes first- order and second- order theories of dissipative fluid dynamics. A 
first-order theory is, e.g., the Navier-Stokes (NS) approximation where the dissipative quantities 
n, q^ (or z/^), and ti^'^ are expressed solely in terms of the primary variables e, p, n, u^, or gradients 
thereof. A second-order theory is, for instance, the so-called Israel- Stew art (IS) theory j9j. Here, 
the dissipative quantities are independent dynamical quantities whose evolution is governed by 
differential equations similar to the fluid dynamical equations ([3]), (j4]), and ([5]). 



2.3 Navier-Stokes Approximation 

In the NS approximation, the dissipative quantities 11, g'^, ■n^^" read: 

Hns = -C^, (6) 

<s = 2r/a'^% (8) 

where (3 = 1/T and a = f3fi; fi is the chemical potential associated with the net charge density n. The 
quantitites (, k, and 77 are the bulk viscosity, thermal conductivity, and shear viscosity coefficients. The 
shear tensor is defined as a^'^ = V^^u'^^. Note that Eq. ([7]) holds in the Eckart frame; in the Landau 
frame, simply replace g^ — z/^(e + p)/n. 

Since 11, q^ (or v^^), and ti^^ are solely given in terms of the primary variables, one can simply insert 
the expressions (EI), ((71), and ([H]) into the Eqs. ([3]), (jll), and ([5]) and obtain a closed set of equations 
of motion: the relativistic generalization of the non-relativistic NS equations. As mentioned in the 
Introduction, the problem with these equations is that they lead to unstable solutions and support 
acausal propagation of perturbations. 



3 



3 Power counting 



In order to derive the IS equations, we first have to discuss the length scales appearing in fluid dynamics. 
Then we identify the Knudsen number as the small quantity in terms of which one can do consistent 
power counting. The IS equations will then emerge at second order in an expansion in powers of the 
Knudsen number. 

3.1 Scales in fluid dynamics 

In principle, there are three length scales in fluid dynamics, two microscopic scales and one macroscopic 
scale. The two microscopic scales are the thermal wavelength, Ath ~ /9, and the mean free path, £mfp ~ 
(((j)n)~^, where (cr) is the average cross section. The macroscopic scale, I/hydro, is the scale over which 
the fluid flelds e, n,u^, . . . vary, i.e., gradients of these flelds are typically of order 9^ ~ -^hydro- Note 
that, since ~ /3 ~ Ath, the thermal wavelength can be interpreted as the interparticle distance. 

Note also that the notion of a mean free path requires the existence of well-defined quasi-particles. In 
strongly coupled theories, the quasi-particle concept is no longer valid. In this case, there are only two 
scales, Ath and Lhydro- 

As a first important result let us note that the ratios of the transport coefficients i^/ Pi and rj to 
the entropy density are solely determined by the ratio of the two microscopic length scales, ^mfp/Ath- 
We demonstrate this explicitly for the shear viscosity to entropy density ratio, similar arguments also 
hold for the other transport coefficients. For the proof, note that r] ~ ((c)Ath)~^ and n ~ ~ s. 



There are two limiting cases, (a) the dilute-gas limit, -^mfp/Ath ~ ^/s — oo, and (b) the ideal- fluid 
limit, £mfp/Ath ~ ?7/s — ^ 0. Estimating ^mfp ~ {a)~^\l-^, the first case corresponds to (cr)/Ath 0, i.e., 
the interaction cross section is much smaller than the area given by the thermal wavelength. In other 
words, the average distance between collisions is much larger than the interparticle distance. In this 
sense, the dilute-gas limit can be interpreted as a weak-coupling limit. Similarly, the ideal-fluid limit 
corresponds to ((j)/A^jj — > oo. This is the somewhat academic case when interactions happen on a scale 
much smaller than the interparticle distance. In this sense, this is the limit of infinite coupling, i.e., the 
interactions are so strong that the fluid assumes locally and instantaneously a state of thermodynamical 
equilibrium. 

For any value ofrj/s (and, analogously, (/s and k/{[3s)) between these two limits, the equations of 
dissipative fluid dynamics may be applied for the description of the system. The situation is particularly 
interesting for £mfp/Ath ~ ?7/s ~ 1 or, equivalently, (cr) ~ A^h ~ T'"^. In this case, there is only a single 
microscopic scale Ath in the problem. This occurs, for instance, in strongly coupled theories without 
well-defined quasi-particles. 

One may derive the equations of dissipative fiuid dynamics in terms of a gradient expansion. How- 
ever, in order to be able to truncate this expansion after a finite number of terms one has to require 
that gradients are small or, equivalently, that the Knudsen number is small. This is discussed in the 
following subsection. 

3.2 Expansion in Terms of the Knudsen Number 

The Knudsen number is defined as K = ^mfp/-f'hydro- Since ^^j^ydro ^ ^ij-i expansion in terms of K is 
equivalent to a gradient expansion, i.e., an expansion in terms of powers of d-raipdn- 

We can now establish a second important result: provided that the dissipative quantities 11, (or 
v^), and 7^^^ do not differ too much from their NS values, the ratios of these quantities to the energy 



Then, 




{cr)n Ath (o-)Ath n s ' 



11 1 1 77 



(9) 
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density are proportional to the Knudsen number. We demonstrate this exphcitly for the bulk viscous 
pressure; similar arguments apply for the heat flux current and the shear stress tensor. We use the 
fundamental relation of thermodynamics, e + p = Ts + nn, to estimate /5 e ~ Ath e ~ s and we employ 
Eq. ([6]) to write 

n IInS C,/) C-^th„ ^„ C [ ^mfp\ ^ u nn\ 

— ~ T Ath 6' -„ — imipOf.u'' — Ku'^r^K. (10) 

e e e Ath e s ^^fp s \\h J 

In the last step, we have employed Eq. ^ and the fact that -u^ ~ 1. The result is remarkable in the 
sense that Il/e is only proportional to and independent of the ratio of viscosity to entropy density. 
The reason is that this ratio drops out on account of Eq. ([9]). 

As a corollary to Eq. ( |T0|) . we conclude that if the Knudsen number is small, K ~ 5 <^ 1, the 
dissipative quantities are small compared to the primary variables. The system is close to local ther- 
modynamical equilibrium or, in other words, close to the ideal-fluid limit. The equations of dissipative 
fluid dynamics can then be systematically and in a well-controlled manner derived in terms of a gradi- 
ent expansion or, equivalently, in terms of a power series in K or, equivalently because of Eq. fllUp . in 
terms of powers of dissipative quantities. At zeroth order in K, one obtains the equations of ideal fluid 
dynamics. At first order in K, one obtains the NS equations. At second order in K, the IS equations 
emerge. 

Finally note that the independence of the ratio of dissipative quantities to primary variables from 
the viscosity to entropy density ratio has important phenomenological consequences. It guarantees 
that, provided that gradients of the macroscopic fluid fields (and, thus, K) are sufficiently small, the 
NS equations are still valid and applicable for the description of systems with large ry/s, e.g. water at 
room temperate and atmospheric pressure. 



4 The Israel- Stewart Equations 

The IS equations of relativistic dissipative fluid dynamics can be derived from the Boltzmann equation 
via Grad's 14- moment method [14j|. The details of this derivation will be presented elsewhere [13j. To 
second order in dissipative quantities (or equivalently, because of Eq. (fTOj) . to second order in K) the 
equations read: 

n = HNs-rnll 

+ Tnqq ■ u - ^iiqd ■ q - C ^^o^O 

+ An, g ■ Va + An. 7r^%, , (11) 

- A,,fT'^'^g, + A,nnV'^a + A,.7r^'^V,a , (12) 
+ 2 T^g q<^'u^> + 2 4g V<'^g^> + 2 ti^^uo^>^ - 2 t] 82 tt'"'' 9 

- 2r.7r;,<V>^-2A.gg<^V">a + 2A.nn(T'^% (13) 
where u'^'^ = \ A''^ (9aM/3 — di^Ua) is the vorticity. We now discuss these equations. 

1. The transport coefficients (, k, t], the relaxation times rn, Tg, r^, the coefficients Tuq, r^n, 'Tqn, 
Tnq, ^uq, ^qU, ^qir, 4g, An,, An,r, A,;,, A^n, Ag., A.g, A.n are (complicated) functions of a, (3, di- 
vided by tensor coefficients of the second moment of the collision integral, for details, see Ref. [13]. 
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As one sends the cross section in the colhsion integral to infinity, all these coefficients go to zero. 
Since then also the NS values vanish, one ends up with the trivial solution U = q^^ = ir^^ = to 
the IS equations, i.e., one recovers the ideal-ffuid limit. This is consistent with the discussion in 
Sec.lSTTl 

2. The coefficients 5o, 5i, 82 are (complicated) functions of a, /3. 

3. The form of the equations is invariant of the calculational frame (Eckart, Landau, . . . ), however, 
the values of the coefficients are frame-dependent. The obvious reason is that the physical inter- 
pretation of the dissipative quantities is frame- dependent. For instance, in the Eckart frame, 

is the heat fiux current, while in the Landau frame, q^ = —v^{e + p)/n is the (negative of the) 
diffusion current, multiplied by the specific enthalpy. Details are given in Ref. [T3] . 

4. The NS terms (the first terms on the r.h.s.) are of first order in K, all other terms are of second 
order in K. Consequently, dropping the latter one obtains the NS equations (jSD, dZl), and ([H]). 

5. The so-called simplified IS equations (in the terminology of Ref. [12]) emerge by keeping only the 
first lines of Eqs. ffTTj) . f|T2|) . and f|T3l) . The resulting equations have the simple interpretation that 
the dissipative quantities 11, q^, and t^^^ relax to their corresponding NS values on time scales 
Tn, Tg, and r^, respectively. 

6. Considering the full IS equations (|TT1) . (fT2l) . and f|T3|) . for times t < Xj, i = 11, g, tt, the dissipative 
quantities 11, g'^, vr'^'^ are driven towards their NS values. Once they are reasonably close to these, 
the first terms on the r.h.s. largely cancel against the l.h.s.. The further evolution, for times 
t > Tj, is then determined by the remaining, second-order terms. These terms therefore constitute 
important corrections for times t > Tj and should not be neglected. 

7. The first two terms in the second line of Eq. (ITT!) , the first five terms in the second line of Eq. 
( IT2|) . and the first three terms in the second line of Eq. (fT3|) were also obtained by Israel and 
Stewart [9] , while the remaining second-order terms were missed or neglected. Presumably, Israel 
and Stewart made the assumption that second-order terms containing 9, a^'^, or V^a are even 
smaller than suggested by power counting in terms of K. The last two terms in the second line 
of Eq. ffTTl) . the last four terms in the second line of Eq. f|T2l) . and the last three terms in the 
second line of Eq. f|T3|) were also obtained by Muronga [llj, while the other second-order terms do 
not appear in that paper. A possible reason is that the corresponding treatment is based on the 
phenomenological approach to derive the IS equations and terms that do not generate entropy are 
absent. The terms in the third line of Eqs. flTT]) . f[T^ . and f[T^ were neither given by Israel and 
Stewart [9j nor by Muronga [11] and are thus genuinely new to this paper (with one exception 
discussed below). 

8. If we set n = = in Eq. (fT3l) . the resulting equation for tt'^'^ is identical to that found in Ref. |16j . 
In particular, the first term in the third line was already obtained in that paper, where it appeared 
in the form [Xi/rj^] 'k^^t:'^^^. Using the NS value ([8]) for vr*^^, which is admissible because we are 
computing to second order in to this order this is identical to 2{Xi/vi) n^'^a^^^. By comparison 
with Eq. (IT^ . we thus get a prediction for the coefficient Ai from kinetic theory, Ai = r^r], in 
agreement with Ref. [16]. Note, however, that this discussion so far neglects additional terms 
which arise at second order in K when expanding the second moment of the collision integral. 
(This was already noted in Ref. [16].) This will change the coefficient of the respective term such 
that it is no longer equal to r^. It will therefore also lead to a different result for Ai. In a recent 
study [T7] a complete calculation was performed. 
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5 Conclusion 



In this paper, we have discussed the full Israel-Stewart (IS) equations of relativistic dissipative fluid 
dynamics as they emerge from applying Grad's 14-moment expansion to the Boltzmann equation and 
truncating dissipative effects at second order in the Knudsen number K = £mfp/-Z^hydro- Our treatment 
is not restricted to the shear stress tensor, but also contains the bulk viscous pressure and the heat flux 
current. It is thus also applicable to non-conformal systems with non- vanishing net charge density. 

We have shown that, in comparison to previous discussions [HI [IHl [H], additional second-order 
terms appear. One of these terms was found already in Ref. [16]. The details of the derivation of the 
full second-order IS equations will be presented elsewhere [13]. The second-order terms are multiplied 
with coefficients whose values depend on the calculational frame. Explicit expressions will be reported 
elsewhere [Hj. Future directions of work comprise the generalization to a system of various particle 
species [15], as well as the numerical implementation and application to modelling the dynamics of 
heavy-ion collisions. 
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